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CO 

'■O We construct a class of anomaly- free supersymmetric U(l)' models that are 

y— I characterized by family non-universal U{1)' charges motivated from Eq embed- 

dings. The family non-universality arises from an interchange of the standard 
roles of the two SU (5) 5* representations within the 27 of E§ for the third gen- 
eration. We analyze U{1)' and electroweak symmetry breaking and present the 
particle mass spectrum. The models, which include additional Higgs multiplets 
and exotic quarks at the TeV scale, result in specific patterns of flavor-changing 
neutral currents in the b —* s transitions that can accommodate the presently 
observed deviations in this sector from the SM predictions. 



1 Introduction 

Extensions of the Standard Model of particle physics (SM) with an additional anomaly- 
free gauged U(l)' symmetry broken at the TeV scale are arguably some of the most 
well-motivated candidates for new physics (for a review, see [I]). Such symmetries 
are theoretically motivated, as they represent the simplest augmentations of the SM 
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gauge sector and are ubiquitous within string and/or grand unified theories. While the 
phenomenology of such Z' gauge bosons depends on the details of the couplings of the 
Z' to the SM fermions, current limits from direct and indirect searches indicate typical 
lower bounds of order 800 — 900 GeV on the Z' mass and an upper bound of ~ 10~ 3 on 
the Z — Z' mixing angle [2]. For a reasonable range of couplings, the presence of such 
TeV scale Z' bosons should be easily discernable at present and forthcoming colliders 
such as the Tevatron and the Large Hadron Collider (LHC). 

Within the context of supersymmetric theories, a plethora of U(l)' models have 
been proposed, including scenarios motivated by grand unified theories (GUTs) such 
as SO (10) and Eg and scenarios motivated from string compactifications of heterotic 
and/or Type II theories (see [1] for a review). Recent models also include scenarios in 
which the U(l)' mediates supersymmetry breaking [BJ, plays a role in the generation 
of neutrino masses [I] and/or spontaneous i?-parity violation [5], or provides a portal 
to a hidden/secluded sector (for reviews, see [HJ Ej). Though the details of the U(l)' 
charge assignments are model-dependent, generically the cancellation of U(l)' anomalies 
requires an enlargement of the matter content to include SM exotics and SM singlets 
with nontrivial U(l)' charges. In these theories, the SM singlets also typically play an 
important role in triggering the low-scale breaking of the U(l)' gauge symmetry. 

In most models of this type, the U(l)' charges of the quarks and leptons are family 
universal. Though this feature is desirable for the first and second generations due to the 
strong constraints from flavor-changing neutral currents (FCNCs), there is still room for 
departures from family universality for the charges of the third generation. In fact, this 
often occurs in string constructions if the families result from different embeddings (see 
e.g., [SI E])- Indeed, though many of the results from the B factories have indicated a 
strong degree of consistency with the Cabibbo-Kobayashi-Maskawa (CKM) predictions 
of the SM, there are hints of non-SM FCNC patterns within the b — > s transitions 
for both AB = 1 and AB = 2 processes at the level of a few standard deviations 
[10J. Of the many options for new physics models that can explain this discrepancy, 
family non-universal U(l)' models are interesting in that they are theoretically well- 
motivated scenarios that lead to tree-level FCNC, as opposed to scenarios in which the 
new physics contributions are loop-suppressed [UJ . A recent model-independent analysis 
of Z'-mediated FCNC in the b — > s transitions showed that this general framework 
can accommodate the data [121 [13]. (Related analyses include [THUS]). However, it is 
optimal to consider the bounds on specific family non-universal U(l)' models in addition 
to the fully model-independent results. 

Our purpose in this paper is to construct and analyze supersymmetric anomaly- 
free family non-universal U(l)' models (which we will denote as NUSSM models). Our 
strategy in building this class of NUSSM models is to exploit the well-known fact that 
in Eq models, there are two options for embedding the down quarks and lepton doublets 
in the 5* representation of SU(5), which is related to the fact that the down-type 
Higgs and the lepton doublets have the same gauge quantum numbers. By choosing 
one embedding for the first and second generations and the alternative embedding for 
the third generation, we can obtain anomaly-free models in which the additional family 
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Table 1: The decomposition of the fundamental 27 representation of Eq with respect 
to the SO (10), SU(5), and U(l)' subgroups. U(l)j and U(l)s correspond to 
specific linear combinations of the U(l) x and U (1)^ gauge groups that result 
from 50(10) -► SU(5) x U(1) X and E 6 -► 50(1) x 17(1)^, respectively. 



£ 6 


5*0(10) 


SU(5) (1st and 2nd families) 


ST/(5) (3rd family) 


2yidQ x 




2Q/ 


2 v / 15Qs 




16 


* 10 = (u L ,d L ,u c L ,e L L ) 


$10 = {tL,b L ,t c L ,T c L ) 


-1 


1 





"1/2 






= (d%v L ,e L ) 


$ 5 . = (A c , 


3 


1 


1 


4 






*i = ^£ 


$! =S 


-5 


1 


-1 


-5 


27 


10 


05 = (A ft u ) 


X 5 = (A,H U ) 


2 


-2 





1 






as. = (-D c , fc,,) 


S 5 * = (b%,v TL ,T L ) 


-2 


-2 


-1 


-7/2 




1 


(T = S 


E = 





4 


1 


5/2 



non-universal C/(l)' is given by a particular linear combination of the usual U(l)^ and 
U(l) x of _E 6 -inspired models. 

This paper is structured as follows. We begin by outlining our basic procedure and 
presenting the resulting classes of anomaly-free family non-universal U(l)' models. In 
the following section, we analyze the gauge symmetry breaking and comment on general 
features of the mass spectrum. We next turn to an analysis of the implications of these 
models for FCNC in the b — > s transitions, then provide our concluding remarks. 

2 £ 6 -Motivated Family Non-universal U(l)' Models (NUSSMs) 

In U(l)' models, the cancellation of gauge anomalies generally implies that additional 
fermions are present in the theory (see e.g., [I]). To motivate the presence of these 
additional fermions and construct simple anomaly-free family non-universal models, our 
approach is to exploit the properties of E e embeddings of the SM fermions and Higgs 
fields in grand unified theories. Recall that in Eq models, the SM particles are embedded 
in the fundamental 27 representations. With respect to the two-step breaking scheme 
of Eq to its SO (10) and SU(5) subgroups 

E 6 -> 50(10) x U(l) f -> SU(5) x U(l) x x U(l)^, (2.1) 

the 27 has the decomposition 

27 = 16 + 10 + 1 = (10 + 5* + 1) + (5 + 5*) + 1, (2.2) 

with respect to the representations of 5O(10) and 577(5), respectively. Hence, the 27 
has two 5* multiplets; these representations are used to embed the down-type SU(2)- 
singlet quarks with the lepton doublets and exotic SU (2)-singlet quarks with down-type 
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Higgs doublets. A standard choice for model-building is to have the down-type quarks 
and lepton doublets of all three SM families in the 5* of the 16, though models with the 
SM down-type quarks and lepton doublets in the other 5* have also been considered in 
the literature [TH [T7] . 

We will assign the down-type quark singlets and lepton doublets of the first and 
second generations to be in the 5* of the 16, and the associated particles of the third 
generation to be in the 5* of the 10, as shown in Table [TJ The matter content of 
these theories thus includes the following fields: (i) the SM first and second families 
{^o,^!*,^}, Higgs plus exotic fields {05,05*}, and singlets a l (i — 1,2 is a family 
index), and (ii) the SM third family {$10, £5*, S }, Higgs and exotics {E 5 ,$ 5 *}, and 
singlet $1. The Higgs sector of the theory thus generically has multiple Higgs doublets 
and sine; lets beyond those of the MSSMQ 

The additional family non-universal U(l)' in these NUSSM models is then a linear 
combination of the U(l) x and the U(l)$ gauge groups: 

Q' = cos 9 Q x + sin 9 (2.3) 

(the assumption is that the orthogonal linear combination of U(l) x and U(l)$ is either 
absent or broken at a high scale). The familiar U(l) v group, which has tan 6 = —a/5/3, 
is family universal and therefore is not useful for our purposes]^] Two viable options for 
the additional U(l)' group are: 

• U(l)i (tan6> = a/3/5). In this model (the inert model) the U(l)' gauge boson 
couplings to the up-type quarks vanish [20] . Hence, the production of the associ- 
ated Z' boson is suppressed at hadron colliders. This is especially the case at the 
Tevatron, since in high-energy pp collisions the Z' production via down quarks is 
suppressed by an order of magnitude relative to up quarks [2T] . 

• U(l)s (tan6> = a/5/27). This symmetry is motivated by models with a secluded 
U(l)' breaking sector and a large supersymmetry breaking A-term that have (1) 
an approximately flat potential that results in an appropriate Z-Z' mass hierar- 
chy [22] ; (2) a strong first order electroweak phase transition and large spontaneous 
CP- violation, which can result in viable electroweak baryogenesis |23j . 

While we use the E e framework to motivate the matter content and U(l)' charges of 
these models, we do not work within a full grand unified theory. More precisely, we do 
not impose the Eq Yukawa coupling relations. This allows for a TeV-scale U(l)' without 
the danger of rapid proton decayj^] 

The allowed superpotential terms of NUSSM models (assuming a conserved i?-parity) 
are the couplings that are consistent with the SM and U(l)' gauge symmetries. An 

MSSM-type gauge unification requires the introduction of an additional non-chiral Higgs pair h+h* 
from an incomplete 27 + 27* [TB] . 

2 This linear combination occurs in certain Calabi-Yau compactifications of heterotic string theory 
if Eg breaks to a rank 5 group via the Hosotani mechanism [19] . 

3 Detailed studies of E e theories with broken Yukawa relations can be found in [TT] . 
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inspection of Table [l] shows that in the language of the SU (5) decomposition, the usual 
10 5*5* and 10 10 5 terms that give rise to quark and lepton Yukawa couplings for 
all three families (including mixing terms between the third family and the other two 
families) are allowed by both U(l) x and U(l)^. Similarly, these symmetries allow the 
generation of Yukawa interactions for the exotic quarks and for the Higgs doublets with 
the SM singlets (i.e., mass terms for the exotic quarks and effective fi terms for the 
Higgs fields, which will be of importance for gauge symmetry breaking). 

For simplicity, we assume that only the neutral Higgs bosons from the third family 
(H ud and S) and one of the first two families (h Ujd and s) acquire vacuum expectation 
values (VEVs)]^] In this limit, the Higgs bosons and Higgsinos in the other family have no 
mixing at leading order with the other particles. The mass eigenvalues of these particles 
are determined by the VEVs of h Uj d, s, H ud , S as well as the Yukawa couplings and soft 
parameters which are not directly involved in the electroweak symmetry breaking. In 
this article, therefore, we will not discuss them in detail. The relevant superpotential 
terms are then given by (I, J = 1, 2, 3 and i,j = 1,2 are family indices): 

W Y = (f^h d + f^H d )Q[di + (h{ J h u + h I 2 J H u )Q I L u J R + 

(fk J h d + fHH^l/ei + (y[ J h u + yl'HJl/vi (2.4) 
W H = X 1 sh d h u + X 2 sh d H u + X 3 SH d h u + X A SH d H u + AW H , (2.5) 

in which the Yukawa couplings satisfy the relations f % d \ = /Jf = 0, f d \ = jffi = 0, 
fH = flf = 0, f% = f* ee 0, yf = yf ee 0, and yf = yf = 0. In Eq. @, AW H 
represents additional superpotential terms that are consistent with U(l)j or U(l)s, but 
explicitly break the orthogonal linear combination of U(l) x and U(l)^. These terms 
are needed to avoid the appearance of undesirable light axions in the low energy theory 
(see [24J for a recent discussion). For the U(l)j model, AWh is a bilinear term: 

AW H = X 5 sS. (2.6) 

Although the coupling A5 in AWh is dimensionful, there is no associated fi problem 
in the traditional sense. This term is not necessary for U(l)' or electroweak symmetry 
breaking, so its mass scale need not be connected with the electroweak scale. The 
Giudice-Masiero mechanism [2S] therefore can be implemented in both gravity- and 
gauge-mediated breaking frameworks to produce such a term, even though the A5 in the 
latter case is typically small. In the U(l)s model, AWh consists of the trilinear term 

AW H = X 5 ssS. (2.7) 

In what follows, we will focus on the U(l)j model as a concrete and minimal example, 
and defer the U(l)s model for future study. 



1 This is actually without loss of generality by appropriate field redefinitions. 
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3 Particle Mass Spectrum and Gauge Symmetry Breaking 



The gauge group of the Z7(l)j NUSSM model is given by SU(3) e x SU(2) L x U(l) Y x 
U(l)i, with gauge couplings g%, g 2 , pi, and (/, respectively. The matter content, which 
was presented in Table [TJ includes three sets of Higgs fields (one pair of doublets and one 
singlet per family) and three sets of exotic down-type quarks in addition to the MSSM 
fields. As previously discussed, we assume that only two of the three Higgs doublet pairs 
(H Ut d and h Ut d) acquire vacuum expectation values, and hence focus on the couplings of 
these Higgs fields only. Restricting to this set of terms, the superpotential for the model 
is given in Eq. Kty, Eq. (pOL and Eq. ggb. 



The tree-level Higgs potential (for the neutral components of the fields) is given by 
V = V F + V D + V H , in which 



\X lS h° u + X 2 sH° u \ 2 + \X lS h° d + X 3 SH° d \ 2 + \X 3 Sh° u + X 4 SH°\ 2 + \X 2 sh° d + X 4 SH° d \ 2 



+\X 1 h° u h d + X 2 H° u h d + X 5 S\ 2 + \X 3 h° u H° d + X 4 H° u H° d + A 5 s| 



(3i 



V, 



G 2 



D 



1 2 



W + \H° U \ 2 - \H°\ 2 ) 2 + 9 - (-K\ 2 + \s\ 2 + \H Q d \ 2 - \S\ 2 ) 



12 



2\ 2 



(3.9) 



+(A Xl X lS h° d h° u + A X2 X 2 sh d H° u + A X[i X 3 SH° d h° u 

+A Xi X 4 SH° d H° u + B Xr X 5 sS + H.C.) , (3.10) 



where G 2 = g\ + g 2 . We also include the one-loop contribution to the potential: 
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1 STYAWm^ 



647T 2 



A^ 
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(3.11) 



in which A4 2 (Hi) denotes the field-dependent mass-squared matrices of the theory, and 
A^g is the MS renormalization scale. We will only consider the dominant one-loop 
contributions that arise from the top quark sector: 
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+ mUHi) In 



A 2 - 

MS 
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(3.12) 



The Higgs potential allows for a rich structure of CP- violating effects, including explicit 
CP violation (for complex couplings) and spontaneous CP violation. In this work, we 
will assume that all couplings are real and let the potential parameters satisfy some 
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necessary constraints such that spontaneous CP violation can be avoided. In this case, 
the vacuum expectation values of the neutral Higgs components can be taken to be real: 

(h° d ) = v 1 , (h° u ) = v 2 , (H° d ) = V u (H° U ) = V 2} (s) = Sl , (S) = s 2 . (3.13) 

Before turning to a numerical analysis, we begin with a general discussion of the particle 
mass spectrum, starting with the gauge bosons. The Z — Z' mass-squared matrix is 

_( Ml M\ z ,\ 
Mz - Z ' ~ 1 Mir, M% I ' (314) 



l zz 



in which 



Ml = ^(vl + vl + V^ + Vi)^^ 



Ml, = 2g'\QZvl + Qlvi + Q'M + QtVi 2 + QtVi 2 + Q'M), 
Ml z , = g'G{Q' hd vl-Q' hu vl + Q' Hd V^-Q' H V 2 \ (3.15) 

with v 2 = v\ + v 1 + V 2 + V 2 2 = (174 GeV) 2 . The mass-squared eigenvalues are 



M kz 2 , 



\ [Ml + M 2 Z , =F v /(Mf-M|,) 2 + 4Ml^ , (3.16) 



and the Z — Z' mixing angle az-z' is 



1 / 2M| Z , 



az - z, = r ictan {Mi,-Mi)' (3 - 17) 

which is bounded to be less than a few times 10 -3 (see [2] for a recent discussion). This 
typically requires that the singlet vacuum expectation values Si i2 3> 1 TeV, resulting in a 
TeV-scale Z' mass. The charged gauge boson mass is given as usual by M w ± = g 2 v/\/2. 
In the basis {B', B, W 3 °, h% h° u , s, H° dl H°, S}, the neutralino mass matrix is 

/M*o(3,3) M,o(3,6)\ 
M *° ~ 1m ? (3,6) t M^o(6,6) ) ' [6AH) 



in which 



m; o o 

M^o(3,3) = | Mi 
M- 



2 



M^o(3,6) 



9iv 2 



V2 



92Vl 



-^ 2 Q2V 2 










hg 2 V 2 



(3.19) 



V2 



7 



M*o(6,6) 



/ A1S1 Aiv 2 + A 2 V 2 A 2 si \ 

A1S1 \ 1 v 1 A 3 s 2 A3V1 

Aiw 2 + A2V2 Ai^i A 2 wi A 5 

A3S2 A4S2 A3W2 + X4V2 

A 2 si A 2 fi A 4 s 2 A4V1 

\ Ag^l A 5 A3W2 + A4V2 A 4 T4 / 



In the above, = y/2g'Q $((])*), and M{, M±, and M 2 are the gaugino mass parameters 
for U(l)', U{\)y, and SU(2) L) respectively. The chargino mass matrix is 

/ M 2 f 2 v 2 f 2 V 2 \ 

M x± = A i fi i A 2 fi i • (3-20) 

A3S2 A4S2 / 

Since 51,2 t>i,2, V\,i because of the experimental bounds on olz-z 1 , the charginos and 
neutralinos are typically heavy unless the A's are small or the gaugino masses are light. 
In the latter situation, the lightest chargino and neutralino will be gaugino-like. 
The mass-squared matrices of the sfermions (denoted collectively as 0) are 

r2_f(Ml)n (M 2 ) 12 



With the definitions 



"i = i$£ $)-.'■ (3 ' 21) 



A ee °^{Tt -QU sin 2 w ){vl - v\ + Vf - V 2 2 ) , (3.22) 

a; ee q;</ 2 (qW + q' h y 2 + q'A + Q< H yt + q' h v? + q> s sD , (3.23) 

the entries for example of the up-type squark mass-squared matrix are: 



(M 2 ) n = ml L + m 2 u + A ilL +A' i 



U L I 



(Mf)ia = /ll(AiUiSi + A3V1S2) + /i 2 (A 2 WlSl + A4V1S2) - (i4 hl /liU2 + A fta /i 2 ^ 2 ) 

(M 2 ) 2 i = (M 2 ) 12 

(M 2 ) 22 = m 2 fl +m 2 +A^ + A' fifl . (3.24) 

Analogous expressions can be written for the down-type squarks, sleptons, and sneutri- 
nos. The physical Higgs spectrum consists of 6 CP-even neutral Higgs bosons, 4 CP-odd 
neutral Higgs bosons, and 6 charged Higgs bosons (not including the second family). 
The tree-level charged Higgs boson mass-squared matrix is given in the Appendix. 

Next we turn to a numerical analysis of this sector of the model, taking into account 
the constraints on the Z' gauge boson. We explore the viable regions of parameter space 
in which (i) s±, s 2 ^> i>i,2, Vi,2, which is needed for a TeV scale Z', and (ii) V 2 > V\ > i>i,2, 
which is motivated by the observed hierarchies in the SM fermion mass spectrum. To 
obtain an acceptable minimum, typically we need the Higgs soft mass parameters to 
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satisfy m 2 s or m| <C m\ u ,m 2 Hu < mj ld ,m 2 Hd . We also set the £7(1)/ gauge coupling to 
|<7i and enforce the following constraints on the Yukawa couplings^ 



hf = hf, hf = hf = hf = hf = 0. (3.25) 
which result in the condition 

, ,„ , oo 165 GeV 

hf = hf = , 3.26 

v 2 + V 2 

in which we have included the one-loop QCD corrections to the top quark mass. 

We consider one typical numerical example; the relevant input parameters and results 
are summarized in Tables [2]-[4j The mass spectrum of the neutral Higgs bosons are 
calculated at one-loop level, and the mass spectra of the other particles are calculated 
at tree- level. As a check, we estimate the Z' mass and the Z — Z' mixing angle by using 
the Higgs VEVs given in Table |2j as follows: 

2 



M Zs « M z , « a/o.18(s? + si) ~ 1.9 TeV, a z . z > « ^fe^ ~ 0.0003, (3.27) 

M z , 

which is consistent with the detailed results. The lightest CP-even Higgs boson H\ is a 
linear combination of the real parts of the four Higgs doublets, with a negligible singlet 
admixture; the orthogonal linear combinations of these four states are the H$, H4, and 
He bosons. These heavier Higgs bosons fall into SU(2) multiplets together with the three 
heaviest CP-odd states and the set of charged Higgs bosons, as follows: (H 3 , A 2 , Hf), 
{Hi, A3, H 2 ), and (Hq, A4, H 3 ). The second lightest and second heaviest CP-even states 
are admixtures of the two singlet Higgs fields, as is the lightest CP-odd boson (which has 
a mass that controlled by the Higgs bilinear terms) . The chargino and neutralino mass 
spectrum is highly model-dependent, as it is sensitive to the electroweak and hypercharge 
gaugino masses, which do not strongly impact the gauge symmetry breaking. Hence, 
the physics of the lightest superparticle (LSP) can vary greatly depending on the exact 
structure of the gaugino sector, though the gauge and Higgs sectors can remain almost 
the same in this case. In our numerical example, in which the gaugino masses are light 
and obey GUT relations, the LSP is a predominantly bino-like neutralino that can be 
an acceptable dark matter candidate in regions of the parameter spacej^] 

Finally, we comment on the exotic colored particles {D % , D a } and {A, A c }. The 
exotic scalars do not obtain VEVs. They and their superpartners influence the gauge 
symmetry breaking only at loop level. These exotic particles are chiral, so their tree- 
level masses can be produced only through Yukawa interactions. The superpotential 



a This approximation must be relaxed slightly to obtain CKM mixing of the third family with the 
first and second families, but that is irrelevant for our present purposes. 

6 The neutralino sector has additional complications due to the presence of the additional Higgs 
supermultiplets that do not participate in electroweak symmetry breaking at tree level. Hence, a 
detailed numerical analysis would be needed to ascertain whether the neutralino LSP satisfies the dark 
matter constraints. As this is tangential to the main purpose of our paper, we do not address it here. 
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Table 2: Parameter values and Higgs VEVs. The dimensional parameter values are 
given in GeV or GeV 2 . The Higgs VEVs are given in GeV. 
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A 4 


A 5 
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0.10 
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0.10 


-0.45 


449 
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Mi 
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m| 
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1.21 x 10 6 
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-1.21 x 10 6 


Am 


Ax 2 


Ax 3 


Ax, 


Bx 5 




-1350 


-1350 


-449 


1080 


-359 


897 


Vl 






V 1 


v 2 


s 2 


54.8 


83.7 


2000 


93.3 


108 


4010 



terms that describe their interactions with the Higgs fields and the corresponding soft 
supersymmetry breaking terms are 

W E = \\ j sD ci D J + X^sETA + ~X] i SA c D i + ~\ 4 SA C A, (3.28) 
V E = A- x fX\ j sD cl D j + A^sD a A + A^\iSA c D l + A~ X4 \±SA c A. (3.29) 

It is straightforward to determine the mass matrices for these states for given Yukawa 
couplings and A parameter values. For 0(Ai ; 2,3,4) ~ 0.1 (where their contributions to 
the effective neutral Higgs potential can be neglected) and not large A x values, the 
exotic particles will typically obtain masses of the order of several hundred GeV. 

4 Z'-mediated FCNC Effects 

In this section, we analyze the Z'-induced FCNC effects After a brief review of the for- 
malism, we will show the results of our correlated analysis of the AB = 1,2 processes 
via b s transitions and discuss the resulting parameter space constraints. The pro- 
cesses of interest include B s — B s mixing and the time-dependent CP asymmetries of 
the penguin-dominated neutral Bd — > (0, V , vr, p,u, fo)Ks decays. 

The FCNC effects in general NUSSM models include both Z'-mediated FCNC pro- 
cesses and contributions to FCNC from the soft supersymmetry breaking parameters. 
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Table 3: The particle mass spectrum and Z — Z 1 mixing angle of the NUSSM model 
(all masses are in GeV). 



M Z2 


sin 9 Z -z> 




m h 


m xt 
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1900 


3.10 x 10" 4 
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m H ± 


633 


1080 
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2340 
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In this work, we assume for simplicity that the soft terms do not result in large FCNC 
effects (this can be easily achieved; see e.g. [26]) and consider only the Z' contributions. 
We now briefly discuss the formalism for addressing such Z' effects in the NUSSM (for 
a model-independent discussion, see [TT 1 H~2l [T3~] ). 

For the SM fermions ipL,R with U{1)' charges fy LR , the fermion mass matrices are 
diagonalized by the biunitary transformation M</, idiag = V^ r M^V^ l (the CKM matrix is 
Vckm = Vu L Vj L ). The chiral Z' couplings in the fermion mass eigenstate basis are 

B*i ee V^Vl , B*» = V4 H e+*V} a . (4.30) 

In our U(l)i model, the only SM fields with nontrivial U(l)' charges are the down-type 
quark singlets and the lepton doublets: 



i LL = ^ I 1 I . (4.31) 





With the unitary matrices V dR ^ L written as 

V dR , LL = f^), (4.32) 

V Y d R ,L L ^d R ,L L J 

where W dR ^ h is a 2 x 2 submatrix, one obtains 

B a R = J_(W\W dR -Y}Y dR W\X dR -Yl R Z dR 

2 \ X d R W d R ~ Z d a Y dn X d R X d R ~ Z d R Z d R 

T i / wt w Lt - y) y Lt wj x L . - y) Z Lj , 

2 \x{W Ll - Z{Y Ll X{X Ll - Z{ l Z Ll > 
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Table 4: The composition of the neutral Higgs mass eigenstates at the one-loop level. 





i 

Kr 


7 n 

Kr 




H° dr 


ttC\ 

ur 


Sf 




0.31 


0.48 


-0.09 


0.53 


0.61 


-0.07 


TT 

H 2 


0.07 


O.Oo 


0.87 


O.Oo 


n nr 

0.05 


0.49 


TT 

^3 


n os 
-U.zo 


U.oo 


n no 
-U.Uz 


n Q/i 
-U.o4 


n oc; 
-U.zo 


n nQ 
U.Uo 


TT 


n as 
-U.oo 


n ii 
-U.ll 


n no 
U.Uz 


n i q 
U. lo 


n /i q 
U.4o 


n n/i 
U.U4 


TT 

-"5 


n nQ 
U.Uo 


-U.U1 


n /in 


n n/i 
U.U4 


n m 
-U.U1 


n Q7 
U.o ( 


TT 


n on 
-U./U 


U.Uo 


n m 
U.U1 


n 7fi 
U. / o 


n fii 
-U.ol 


n nQ 
-U.Uo 




i 


i o 


Si 


Hi 




Si 




-0.04 


-0.02 


0.89 


-0.01 


-0.01 


0.45 


A 

A 2 


0.26 


0.87 


0.03 


0.35 


-0.24 


0.02 


A 3 


0.89 


-0.10 


0.04 


-0.12 


0.43 


0.01 


A, 


-0.20 


-0.08 


-0.01 


0.76 


0.62 


0.02 


Gi 


0.30 


-0.33 


0.20 


0.59 


-0.52 


-0.39 


G 2 


-0.12 


0.21 


-0.40 


-0.25 


0.26 


0.81 



Table 5: The composition of the charged Higgs mass eigenstates at tree level. 





K 


h + * 


H d 




HT 


0.26 


0.87 


0.34 


-0.23 


H 2 ~ 


0.89 


-0.10 


-0.12 


0.42 




-0.20 


-0.07 


0.75 


0.62 


Gx 


0.31 


-0.47 


0.55 


-0.62 
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To avoid the constraints on non-universality for the first two families from K — K mixing 
and fi — e conversion in muonic atoms, we assume small fermion mixing angles or small 
Xd R ,L L , Yd R ,L L elements. The Z' couplings then take the form 

B£,B%n± B&,B%~0{X Ll ,Y Ll ). (4.34) 

Here Bfg' LL and B dR ' Ll generically are complex. The Z'-induced corrections to the 
Wilson coefficients in the U(l)j modeQare given by (for the associated operators, see 

e.g. my- 

AC*f s = -(B^) 2 , AC 3 = - TT^ BfeB^, AC 9 = 7^ B^ s Bf d , 

AC W = -AC 10A = --^—B^. (4.35) 

To achieve sufficient precision, we need to have an accurate knowledge of the relevant 
Wilson coefficients at the b quark mass scale m;, = 4.2 GeV (for general discussions, see 



e.g. [27J). The parameter values used in our calculations are summarized in Appendix B 



• B s — B s mixing. The new physics (NP) contributions to the off-diagonal mixing 
matrix element are parametrized as 

Af£ = (M*°) SM C Bs e 2i <, (4.36) 

where Cb s = 1 and 0^ = in the SM limit. Although the data indicate that Cb s — 1, a 
recent analysis [10] suggests that deviates from zero at the 2 — 3a level (see Table (6j); 
an earlier discussion was given in [28]. The analysis of [10] includes all available results 
on B s mixing, including the tagged analyses of B s — > ip<p by CDF [2H] and D0 [30J. As 
discussed for example in [3T], this discrepancy disfavors scenarios with minimal flavor 
violation (MFV), though no single measurement yet has a 3a significance. 
In our U(l)j NUSSM model, Cb s and 0^" at the m& scale are given by 

C Bs e 2i4B ° = 1 - 3.59 x 10 5 (ACf s + ACf s ) + 2.04 x 10 6 A(5f s . (4.37) 



The large coefficients of the correction terms in Eq. (4.37) are due to the fact that the 



NP is introduced at tree- level while the SM limit is a loop-level effect. 



• B d — > (-0, 7r, <fi, r)',p, ou, f°)Kg decays. The direct and the mixing-induced CP asym- 

7 In the U(l)' s model, in which all SM fermions are charged under the U(l)' symmetry, the Z'- 
induced corrections to the Wilson coefficients take a more general form (see e.g. [321 03]). We will not 
discuss them here. 
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Observable la C.L. 2a C.L, 



m°] (si) 

Ctl (S2) 



-20.3 ± 5.3 
-68.0 ± 4.8 
1.00 ± 0.20 



[-30.5,-9.9] 
[-77.8,-58.2] 
[0.68,1.51] 



Table 6: The fit results for the B s —B s mixing parameters [TO] . The two solutions 
("SI" and "S2") result from measurement ambiguities; see [TO] for details. 



metries in hadronic Bd decays are parametrized as follows: 

1-|A/ CP | 2 <, 2Im[; 

l + |A /cp r *°* l + |A /o 

in which 



P I" I A/ C p | 2 21m [Afcj 

C /cp - TTT\ vi > 6 /cp - T~TT\ x%- \^.iV>) 



*fc P =Vfc P e- 2l<t>B ^. (4-39) 

A f C P 

Here 0e d is the Bd — Bd mixing angle, Af CP is the decay amplitude of Bd — > fcp (^4/ CP 
is its CP conjugate), and r]f CP = ±1 is the CP eigenvalue for the final state fcp- In the 
SM, (p Bd — (3 = arg [— (V c dV^) / (y t dV t l)}, and a non-trivial weak phase enters Af CP only 
at 0(X 2 ). This implies the following SM relation between the decays proceeding via 
b — > sqq(q = u,d,c,s) and the penguin-dominated modes such as Bd — > (n,(j),r}' ,tt, p, 

~ Vf CP Sf CP = sm 2/5 + 0(A 2 ), C /cp = + 0(\ 2 ). (4.40) 

However, the experimental values of sin 2/3 obtained from the penguin-dominated modes 
are below the SM prediction and the results from the charmed Bd — > ipKs mode. The 
central values of the direct CP asymmetries of Bd — > <pKs and Bd —* ujK$ are also 
small compared to the B d — > mode, as shown in Table [7j Since B d — > ^K s is 

a tree-level process in the SM, large values for ASf CP = —i]f CP Sf CP + rj^KsS^Ks an d 
^Cfc P = Cfcp ~ C-ipK s may indicate the presence of NP in the b — > s transitions. 

In NUSSM models, Z'-induced FCNC effects can provide dramatic changes to the 
results, since a new weak phase can enter Af CP at tree level. Following [33], the \f CP 
parameters of B d — > (ip, <f>, rf, ir, p, u, f°)K s at the scale are given by 

X^k s = (-0.63 + 0.74z) (4.41) 
[1 - (2.93 - 2.61i)(AC 3 + AC 3 )* - (2.94 - 2.95z)(AC 5 + AC 5 )* 
+(0.18 - 0.0H)(AC 7 + AC 7 )* - (0.06 - 0.04i)(AC 9 + AC 9 )*] 
/[I - (2.80 + 2.61z)(AC 3 + AC 3 ) - (2.74 + 2.99i)(AC 5 + AC 5 ) 
+(0.17 + 0.01z)(AC 7 + AC 7 ) - (0.04 + 0.05i)(AC 9 + AC 9 )\ , 



14 



fcp 


-VcpSfcp C.L.) 


C fap (la C.L.) 


ijKs 


+0.672 ± 0.024 


+0.005 + 0.019 


cf>K s 


,n 44+0.17 


-0.23 + 0.15 


rjKs 


+0.59 ±0.07 


-0.05 + 0.05 


7CK S 


+0.57 + 0.17 


+0.01 + 0.10 


pK s 




-0.01 + 0.20 


loK s 


+0.45 ± 0.24 


-0.32 + 0.17 


f K s 


+0.62i°;^ 


0.10 + 0.13 



Table 7: The world averages of the experimental results for the CP asymmetries in 
Bd decays via b — ► qqs transitions [32] . 

\ nKs = (-0.70 + 0.70*) (4.42) 
[1 - (1.09 + 0.50*)(AC 3 + AC 3 )* - (6.73 + 2.79i)(AC 5 + AC 5 )* 
-(9.68 + 3.21?)(AC 7 + AC 7 )* + (13.86 + 4.48i)(AC 9 + AG,)*] 
/[I + (1.08 + 0.48z)(AC 3 + AC 3 ) - (6.66 + 2.70«)(AC 5 + AC 5 ) 
-(9.58 + 3.09z)(AC 7 + AC 7 ) + (13.71 + 4.31i)(AC 9 + AG,)] , 



A^ s = (_o.70 + 0.70i) (4.43) 
[1 - (28.62 + 11.37i)(AC 3 + AC 3 )* - (24.08 + 10.4M)(AC 5 + AC 5 )* 
+ (14.57 + 5.88i)(AC 7 + AC 7 )* + (15.08 + 5.92z)(AC 9 + AC 9 )*} 
/[l - (28.27 + 10.89i)(AC 3 + AC 3 ) - (23.80 + 10.00i)(AC 5 + AC 5 ) 
+ (14.39 + 5.64i)(AC 7 + AC 7 ) + (14.90 + 5.67i)(AC 9 + AC 9 )} , 



X n 'K s = (-0.70 + 0.69i) (4.44) 
[1 - (10.88 + 3.29«)(AC 3 + AC 3 )* + (8.26 + 2.06z)(AC 5 + AC 5 )* 
+(2.11 + 0.67z)(AC 7 + AC 7 )* + (2.10 + 0.54z)(AC 9 + AG})*] 
/[l - (10.73 + 3.21i)(AG 3 + AC 3 ) + (8.14 + 2.00z)(AC 5 + AG 5 ) 
+(2.08 + 0.65i)(AC 7 + AC 7 ) + (2.07 + 0.52i)(AC 9 + AG 9 )] , 



\ pKs = (-0.74 + 0.65i) (4.45) 
[1 + (0.26 + 0.062)(AC 3 + AC 3 )* - (19.62 + 1.8H)(AC 5 + AC 5 )* 
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'5, 



-(39.11 + 3.31z)(AC 7 + AC 7 )* - (48.28 + 4.12z)(AC 9 + AC 9 )*] 
/[l + (0.25 + 0Mt)(AC 3 + AC 3 ) - (19.28 + 2.792)(AC 5 + AC 5 ) 
-(38.46 + 5.282) (AC 7 + AC 7 ) - (47.48 + 6.55i)(AC 9 + ACg)] , 

Kk s = (-0.71 + 0.70?) (4.46) 
[1 + (90.48 + 13.54i)(AC 3 + AC 3 )* + (85.24 + 12.50i)(AC 5 + AC 5 )* 
+(32.21 + 4.80i)(AC 7 + AC 7 )* + (19.07 + 2.79i)(AC 9 + ACg)*] 
/[l + (90.01 + 13.29i)(AC 3 + AC 3 ) + (84.80 + 12.26i)(AC 5 + AC 5 ) 
+(32.04 + 4.7H)(AC 7 + AC 7 ) + (18.97 + 2.74z)(AC 9 + ACg)}, 

\ f0Ks = (_o.70 + 0.70i) (4.47) 
[1 + (1.02 + 0.42z)(AC 3 + AC 3 )* - (1.67 + 0.972)(AC 5 + AC 5j 
+(3.19 + 0.93z)(AC 7 + AC 7 )* - (0.12 + 0.15z)(AC 9 + ACg) 
/[l + (1.01 + 0.40z)(AC 3 + AC 3 ) - (1.65 + 0.95*)(AC 5 + AC 5 
+(3.16 + 0.90z)(AC 7 + AC 7 ) - (0.12 + 0.15i)(AC g + ACg)} . 

These results are more general than those of [121 [13], as they include the Z' contributions 
to both the QCD and electroweak penguins. At the leading order, the deviations for 
Cf CP and Sf CP from their SM predictions are a linear combination of these two classes 
of Z' contributions. This discussion is independent of the details of the U{1)' charges, 
so it can be applied to other family non-universal models as well; however, in the U(l)i 
model, the only non-trivial corrections are AC 3 and ACg. 

We now turn to a numerical analysis of the FCNC constraints with the ^(1)/ model, 
for which the three free parameters are \B^\, <p^ s and B^ d . First, we consider B s — 
B s mixing, which involves two of these parameters, \B^ S \ and (f)f s . The experimental 
constraints on these two parameters are illustrated in Fig. [TJ where the various colors 
of the points specify the different confidence levels (C.L.) that the relevant Cb s and 
values represent. There are two separate shaded regions in this figure. The left 
one corresponds to the solution "SI" and the right one corresponds to "S2" (see 
Tab. [HJ. (f)f s varies within the ranges —80° ~ —20° and —90° ~ —70° in the two regions, 
respectively. This is similar to what happens to <p^ s in the LL limit in [T3] , since the AC^ S 
contributions to CB s e 2t ^ Bs in Eq. (4.37) are absent in both cases. In addition, to explain 
the observed discrepancy in B s — B s mixing from the SM prediction, \B^ S \ is required 
to be ~ 10~ 3 . As discussed in [121 [13], there are two reasons for this feature. First, Cb s 
does not deviate significantly from its SM prediction (the anomaly in B s — B s mixing 
is mainly caused by the phase <f>g P )- Second, the corrections of a family non-universal 
Z' arise at tree level, so only a small coupling is needed to explain this small deviation, 



according to Eq. (4.37). The smallness of \B^ S \ is consistent with our assumption of 
small fermion mixing angles, since B^ s is proportional to them (see Eq. (4.34)) as well 
as to g'Mz 1 /{gzMz 2 ). The constraints from the branching ratio Bt(B s — > /i + /i~) and 
Br(5d — > n + n~) can be easily satisfied due to the smallness of |-B^| |13j . 
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Figure 1: Correlated constraints on \B^\ and (j)^ s . Random values for Cb s and 

from the experimentally allowed regions at different C.L. (see Table [TJ are 



mapped to the \B^ S \ 



plane using Eq. (|4.37|). 



With the constrained values of \B£,\ and 4>^ s by B s — B s mixing, we illustrate the 
NP contributions to C^i^^uj^Ks anc ^ ^{<t>,ri'^,p,^Jo)K s m Fig. |2j In this case, the third 
parameter (B dd ) is also involved. For the channel Bd —>■ nKs, we take a strategy 
different from that used in [121 IE], in which it was assumed that the NP enters the 
hadronic decays of neutral Bd meson only through electroweak penguins. In that case, 
the NP effects in the Bd — > irKs channel can be resolved into a factor qe 1 ^ [M]; the 
constraints on this factor from a \ 2 fit °f B — > 7rifg and B — > tttt data have been 
studied in [35J. For our NUSSM model, the NP enters generically through QCD as 
well as electroweak penguins, and hence we treat this channel in the same way as the 
other Bd decay channels. We also assume a 15% uncertainty in the SM calculations 
for each of these modes and a 25% uncertainty for the NP contributions. Here 15% is 
a typical uncertainty level for the hadronic matrix elements of the SM FC operators 
(see e.g. [36]) that is needed to explain the experimental results for C^k s and S^k s 
in the SM [121 C3]. The difference of the uncertainty levels between the SM and 
NP calculations arises because the hadronic matrix elements of the FC operators in 
the SM are better understood than those of the NP operators. To see whether the 
anomalies in B s — B s mixing and the Bd — > ((f), T]',7r,p, u,fo)Ks CP asymmetries can 
be simultaneously accommodated, we have carried out a correlated analysis within the 
U(l)j model. The distributions of \B^\, and B dd constrained at different C.L. are 
illustrated in Fig. [3} Indeed, there exist parameter regions for which the tension between 
the observations and the SM predictions are greatly relaxed. 

In Figs.§and|3| we require 0.005 < \B^\ < 0.5. Given that \B% d \ « g' M Zl /(2g z M Z2 ), 
we immediately find that 1 TeV < Mz 2 < 10 TeV for g' ~ gz- Here B dd can be positive 
or negative, since it resolves a minus sign from the degeneracy of two solutions in B^ s 
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that is specified by a 7r phase difference [T2| H~3] . The red lines in Fig. [3] represent the 
parameter region discussed in our numerical example in which \Bf d \ ~ 0.02. Indeed, we 
see that the anomalies in the hadronic B^ meson decays can be explained simultaneously, 
given the B^ s values required to fit the B s — > B s mixing data. 

5 Discussion and Conclusions 

In this paper, we have discussed a class of family non-universal U(l)' models based 
on non-standard E§ embeddings of the SM that interchange the standard roles of the 
two 5* representations present in the fundamental 27 representation of Eg for the third 
family. The NUSSM models in this class are simple and anomaly-free. They are not full 
E 6 grand unified theories, so the £7(1)' breaking can occur at the TeV scale, resulting 
in a TeV-scale Z' gauge boson that can mediate FCNC in the b — > s transitions. We 
analyzed a representative example of a NUSSM model (the U(l)i model), in which we 
described the low energy spectrum of the theory and determined the constraints on the 
family non-universal Z 1 couplings from the B sector. NUSSM models such as the U(l)i 
model are characterized by a rich spectrum of states with masses at the electroweak to 
TeV scale. The Z'-mediated FCNC in the U(l)i model can easily accommodate the 
observed discrepancies in the b — > s transitions. Related observables such as r — > \i and 
r — > e can also be studied in NUSSM models; we defer this to future work. 
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Appendix A Tree-level Mass-squared Matrix for Charged Higgs 
Bosons 

For charged Higgs bosons, the entries of its mass-squared matrix M^ ± at tree level are 
given in the basis {h d , h+\ H~,H+*} by 



(M 2 H ±)n 


= -^(H 2 -KI 2 + N 2 -|^r) + f(N 2 + N 2 


-hi 2 ) 




(-KI 2 + \si\ 2 + l^il 2 - N 2 ) + (h| 2 + |A 2 | 2 )| Sl | 2 + m 2 hd 


{M 2 H± ) 22 


= % (H 2 - hi 2 + N 2 - |vi| 2 ) + f (hi 2 + hi 2 - 

+ |A 1 | 2 |s 1 | 2 + |A 3 | 2 |s 2 | 2 + m 2 ii , 


• hi 2 ) 




= ~ (hi 2 - hi 2 + \V 2 \ 2 ~ N 2 ) + 9 pv 2 \ 2 + \V 2 \ 2 


-hi 2 ) 




+ Y (-hi 2 + N 2 + l^il 2 - hi 2 ) + (N 2 + N 2 ): 


s 2 | 2 + m 2 Hd 


(M| ± ) 44 


= ^ (hi 2 - hi 2 + \v 2 \ 2 - \v,\ 2 ) + |(M 2 + \v,\ 2 - 

+ |A 2 | 2 | Sl | 2 + |A 4 | 2 |s 2 | 2 + m^ , 


• hi 2 ) 


(M H± )i2 


= (M 2 H± )* 21 = 9 \v\v\ - X x {X\v\v\ + \*V*vl + \;s* 2 ) - 


AiAiSi , 


« ± ) 13 


= (M 2 H± ); 1 = 9 Xiv 1 + (\ 1 x; + x 2 xi)s 1 s* 2 , 




(M 2 H± ) U 


= (M 2 H± )li = ^vlV* - A 2 (A>X + X*V*vl + Xls* 2 ) - 


-4a 2 A 2 si , 


{M 2 H± ) 23 


= (M 2 H ±)* 32 = |Wi - A*(A 3 /^ + X 4 V 2 V X + A 5 si) - 


A* 


(M 2 ±) 24 
(M 2 ± ) 34 


= (M 2 H± )l 2 = ^v 2 V 2 * + A*A 2 si 2 + A*A 4 |s 2 | 2 , 

= (m 2 ± ): 3 = 9 h/*v* - a 4 (a^ 2 x + x\v;v* + a^) 


- ^A 4 ^4S 2 



These entries can be applied to both cases with and without CP violation. 
Appendix B Parameters 

The parameters used in our numerical analysis are summarized below: 



(1) QCD and EW Parameters 

G F = 1.16639 x 1(T 5 GeV- 2 , = 225 MeV, 

M w = 80.42 GeV, sin 2 9 W = 0.23, 
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T)2B = 0.55, 
a s (M z ) = 0.118, 
A = 0.2252, 

p = 0.145, 

Rb = y fp^f 



J 5 = 1.627, 

a em = 1/128, 
A = 0.8117, 
fj = 0.339, 
0.378. 



(2) Masses, Decay Constants, Hadronic Form Factors and Lifetimes 



M w ± = 0.139 GeV, 
M K = 0.498 GeV, 
M = 1.02 GeV, 
Mrf = 0.958 GeV, 
M p = 0.776 GeV, 

Mfo -- 
X v = 0.57, 
m u (pL = 4.2 GeV) = 1.86 MeV, 
m s {n = 4.2 GeV) = 80 MeV, 
m b (fi = 4.2 GeV) = 4.2 GeV, 
h = 237 MeV, 
U = 130 MeV, 
U = 410 MeV, 
f p = 209 MeV, 
F^(0) = 0.330, 
Ff K (0) = 0.379, 
F B/ (0) = 0.250, 

= 0.280, 
t b q = 1.530 ps, 
M Bs = 5.37 GeV, 



M n o = 0.135 GeV, 
M B = 5.279 GeV, 
My, = 2.097 GeV, 
M w = 0.783 GeV, 
M„ = 0.548 GeV 



0.980 GeV, 



Y v = 0.82, 

m d (n = 4.2 GeV) = 4.22 MeV, 
m c (/i = 4.2 GeV) = 0.901 GeV, 
m t (fi = M z ) = 171.7 GeV, 
f B = 190 MeV, 
f K = 160 MeV, 
f u = 200 MeV, 
ffo = 180 MeV, 
Ff K (0) = 0.391, 
Af w (0) = 0.280, 
f/*(0) = 0.030, 



f Bs \jB Bs = 0.262 
tb- = 1.65 ps, 
r Bs = 1.47 ps, 
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Figure 2: The NP contributions to C^,^ ) p )U)> fo ) Ks and S^ jV >^ )P)Wj f ^ Ks , with \B^ S \, 
(f)P constrained by B s — B s mixing. The colors specify the C.L. that their 
inverse image points represent in Fig. [T] (yellow for lex and blue for 2a). 
The boxes specify the allowed regions at la and 1.5a, and the dark points 
denote the SM limit. 
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Figure 3: The \B, ' \, 4> b ^ and B^ d distributions, with values constrained by B s — B s 
mixing at xa C.L. and selected by C ( ^/ ^ PtU j )K s and S (M ■ ,-k, p ,uJ )k s at 
ya C.L.. Here x = 2.0 and y = 1.7 for the purple points, x = 1.0 and 
y = 1.7 for the blue points, x = 1.0 and y = 1.4 for the dark points. The 
red lines represent the vacuum considered in Section [3] in which the values 
of \B^ S \ and (f)^ s are not fixed. 
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